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Abstract. Motivated by the first detections of rotation signatures in the DG Tau jet (Bacciotti et al. 2002 1, we examine possible 
biases affecting the relation between detected rotation signatures and true azimuthal velocity for self-similar MHD disc winds, 
taking into account projection, convolution as well as excitation gradients effects. We find that computed velocity shifts are 
systematically smaller than the true underlying rotation curve. When outer slower streamlines dominate the emission, we 
predict observed shifts increasing with transverse distance to the jet axis, opposite to the true rotation profile. Determination 
of the full transverse rotation profile thus requires high angular resolution observations (< 5 AU) on an object with dominant 
inner faster streamlines. Comparison of our predictions with HST/STIS observations of DG Tau clearly shows that self-similar, 
warm MHD disc wind models with A = 13 and an outer radius of the disc ^ 3 AU are able to reproduce detected velocity shifts, 
while cold disc wind models (A > 50) are ruled out for the medium- velocity component in the DG Tau jet. 
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1. Introduction 

It is now widely accepted that a large scale magnetic field is 
responsible for both the acceleration and collimation of jets 
around young accreting stars. However, the exact launching 
zone (stellar surface, disc truncation radius or wide range of 
disc radii) remains subject to debate. So far, studies aimed 
at constraining proposed ejection models have concentrated 
on the jet collimation, poloidal velocities, and excitation con- 
ditions (Dougados et al. I2000I Lavalley-Fouquet et al. [2000 
Bacciotti et al. l200Q| Woitas et al. l2002l . 

The detection of rotation signatures in the DG Tau jet 
has recently opened new prospects to constrain MHD ejec- 
tion models (Bacciotti et al. 2002). In particular, approximate 
launching radii 0.3-3 AU were inferred in DG Tau, indica- 
tive of a disc wind (Anderson et al. 2003 1. In this Letter, we 
investigate two classes of self-similar disc wind models (Sect. 
2) and analyse in detail the biases introduced in their observed 
rotation signatures by projection, beam dilution and emissivity 
gradients within the jet (Sect. 3). Implications for a proper in- 
terpretation of the DG Tau observations are discussed in Sect. 
4. 
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2. Disc wind models 

We concentrate our study on two classes of steady, self-similar 
MHD keplerian accretion discs driving jets. They are mainly 
characterized by 3 nondimensional free parameters (see Cabrit 
et al. 11999ft : e = Ii/Rq ^0.1, where h is the disc scale height 
at the disc radius Rq, a m - 1, related to the magnetic diffu- 
sivity parameter, and £ = d(logM acc )/d(log/?), which controls 
the mass loading onto field lines. In these solutions, £ is related 
to the magnetic lever arm A (R a /Rq) 2 (R a is the cylindri- 
cal radius at the Alfven surface) by the relation A = 1 + 1 / (2£) 
(Casse & Ferreira 2000a) and therefore A does not vary with Ro. 
Ferreira ( 1997 ) obtained cold solutions powered by the Lorentz 
force only, i.e. where enthalpy is negligible (resulting in typical 
A values 50). Detailed comparison with recent observations 
showed that these solutions reproduce the collimation proper- 
ties of T Tauri microjets (Dougados et al. 2000), but give ex- 
cessive terminal poloidal velocities (Garcia et al. 2001 :>). More 
recently, Casse & Ferreira (2000b) computed warm disc wind 
solutions where entropy is injected at the flow base allowing 
for larger mass loads (A 10) hence lower asymptotic veloci- 
ties. In the following, we adopt a cold solution with A = 50 and 
a warm solution with A = 13. Dimensional scaling parameters 
are set as follows: the inner and outer radii of the disc involved 
in the ejection process are R; = 0.07 AU (typical corotation 
radius for a T Tauri star) and R e = 1 AU (where molecules are 
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Fig. 1. Velocity components (V P ,V^) in the (^,Z) plane for the 
warm solution {A = 13). Dashed lines: varying from 4 to 
16 km s _1 by steps of 4 km s _1 . Greyscale: V p . 



supposed to start to form), the accretion rate through the disc 
M acc = 10~ 6 M yr~' and the central stellar mass M* = O.5M . 
The ejection efficiency (2Mj et /M acc = £ x lnR e /Ri) is then 3 % 
for the cold solution and 10 % for the warm one. 

Asymptotic toroidal velocities in MHD disc winds depend 
mainly on the magnetic lever arm A and Rq. In Fig. [2 we show 
poloidal and toroidal velocities in the (R,Z) plane for the warm 
solution. The cold solution shows a similar map with predicted 
toroidal and poloidal velocities larger by a factor 2. For a given 
streamline, rapidly decreases with distance Z along the jet 
until the maximum expansion radius is reached (at Z/R - 25 
and 20 for the cold and warm solutions respectively). 

To calculate predicted rotation signatures for the above 
models, we constructed synthetic long-slit spectra perpendicu- 
lar to the jet axis (transverse PV diagrams) and convolved them 
by a two-dimensional Gaussian beam to simulate the instru- 
mental spatial and spectral resolutions. The effect of rotation 
is to induce a "tilt" in transverse PV diagrams (dashed lines, 
Fig. Ek). We then get a synthetic velocity shift V^at by cross- 
correlating extracted velocity spectra at symmetric positions 
with respect to the jet axis (Fig. c). Our computations show 
that V s hift values are only weakly dependent on the adopted 
spectral resolution, provided line profiles are well sampled. We 
will use here the a velocity resolution of 50 kms _I , which is 
accessible to current imaging spectrographs, and smaller than 
line profile widths. We adopt a typical inclination of i = 45° 
with respect to the line of sight (Pyo et al. 2003 1 and consider 
transverse PV diagrams at Z pro j = 50-60 AU. 

We will use the [O i] /16300 line, the effect of a tracer with 
lower critical density being discussed later. We explore spatial 
resolution effects by varying the beam size (FWHM) between 
1 AU and 14 AU. To calculate line emissivities in the cold so- 
lution, we use the thermal and ionization structure computed 
a posteriori by Garcia et al. (2001 1) with ambipolar diffusion 
heating. The resulting T e is roug hly constant (T e ~ 10 4 K) 
and the ionization stratification is close to a x e oc \/R law. 
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Fig. 2. a) Synthetic transverse Position- Velocity (PV) diagram 
for the warm solution (integrated over Z pl0 j = 50-60 AU) in 
the [O i] /16300 line, convolved with a 14 AU X 50 kms" 1 
beam, i = 45° and x e = 0.1 x (O.lAU/^o). Dashed lines are 
indicative of peak velocity shifts at 10 and 20 AU from the 
jet axis, b) Velocity profiles extracted from PV diagrams at 
R plo j = ±10 AU, convolved with a spatial beam of 14 AU for 
x e = 0.1 x (0.1AU/#o) (solid lines) and x e = 0.1 x (# /lAU) 
(dashed lines), c) and at R pro j = ±20 AU, with spatial beams of 
1 (dotted) and 14 AU (solid) for x e =0.1x (0.1 AU/R ). 



This gives line profiles with a strong high velocity component 
(HVC), originating from the inner faster streamlines (Garcia et 
al. l2001b V In the case of the warm solution, where the full 
thermal solution is not yet available, we keep T e = 10 4 K 
and we explore two extreme ionization fraction laws: one sim- 
ilar to the cold solution, x e = 0.1 x (0.1AU//?oX where the 
HVC dominates (solid curves in Fig. |2J>), and one with x e = 
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Fig. 3. Velocity shifts predicted by the disc wind models as a function of projected transverse radius /? plo j. Left: cold solution 
with the thermal structure of Garcia et al. (2001 1). Right: warm solution with two ionization fraction laws (see Sect. 2 for 
more details). V s hif t are obtained by cross-correlating [O i] /16300 line profiles at +^ p ,oj, extracted from a transverse PV diagram 
at Z pro j = 50-60 AU. Solid and dash-dotted lines: V s hif t computed using the entire velocity profile, short-dashed and long- 
dashed lines: or restricting to the MVC. Dotted lines show the projected underlying true azimuthal profile extracted from the 
corresponding wind solution. Rj et is the outermost jet radius (for R e - 1 AU). 



0.1 x (Rq/IAU), where emission from outer slower streamlines 
is enhanced, and the profile is dominated by a medium veloc- 
ity component (MVC, dashed curves, Fig.f2J>). The latter case 
reproduces better observed line profiles in the DG Tau jet (see 
Bacciotti et al. 2002 1, and allows us to investigate the effect of 
various emissivity gradients in the jet on the observed V S Mft- 

3. Detailed analysis of disc wind predictions 

We first discuss V s hift values obtained by cross-correlation over 
the entire velocity profile, plotted in Fig. We also plot for 
comparison as a dotted curve the projected azimuthal veloc- 
ity profile ( V'^ftCK) = 2 x V^R) x sin i) extracted from the 
corresponding wind solution at the same Z pro j. We find that 
measured velocity shifts always underestimate the true rota- 
tion profile, especially at small transverse distances from the jet 
axis. This underestimate comes from the integration of the flow 
along the line of sight, combined with projection effects: At a 
given radius R pio j from the jet axis, the line of sight intersects a 
range of flow surfaces with true cylindrical radii R > R pl0 y Due 
to the keplerian law, they rotate more slowly than V,f,(R pr0 j). 
In addition, these outer surfaces are not tangent to the line of 
sight and their projected rotation speed is sin i x V$(R) x sin <p 
with sin <f> = R pro j/R introducing an additional reduction factor. 
However, the exact amount by which the rotation profile is un- 
derestimated by these projection effects critically depends on 
excitation gradients and beam dilution, as we now discuss. 

We first consider the case where inner streamlines (HVC) 
dominate line profiles, which is illustrated by the cold solution 
and the warm solution with x e oc \/R (solid curves, Fig.fJJ. A 



given beam tends to decrease V s un at projected radii below the 
beam diameter (see 14 AU beam, Fig.|5J because of cancella- 
tion effects between opposite sides of the jet. When R pro j varies 
from the beam diameter to the jet outer radius Rj et , each profile 
is dominated by the surface at R — R pro j and the discrepancy 
from the theoretical value is small. Only at radii smaller than 
the central hole radius Z?hoie (4 AU for both solutions), V s hift 
tend to zero as R pro j/Rhoie, due to the above projection effect. 

We examine now the case where outer streamlines (MVC) 
dominate line profiles, which is illustrated by the warm solu- 
tion with x e oc R (dash-dotted curves, Fig. [3}- Line profiles at 
R plo j are now strongly contaminated by surfaces with true radii 
R > R pr0 j, so that y s hift lies significantly below V$S (a factor 
2 at 10 AU), regardless of beam size. Only V s hift measured at 
projected distances on the order of Rj et fit correctly VijJJ, the 
projection factor R pr0 j/R being then ~ 1. 

Finally, when cross correlation is made on the MVC only 
(long-dashed and short-dashed curves, Fig. [3j, Vshift is even 
lower and varies almost as /?prqj /-Rjet- This can be understood 
by considering that cross-correlation on the MVC is sensitive 
only to a narrow range of outer streamlines at R ~ Rj et . Beam 
smearing and excitation gradients have little effect because this 
region is far from the jet axis and has roughly homogeneous 
excitation conditions. 

We have verified that the above results remain valid over a 
broad region of the parameter space. Variation of M*, Ri and R e 
only introduces scaling factors (V oc sJM^/Ro), R^ \ e and Rj et 
being about proportional to R[ and R e respectively. V s hift ob- 
tained from a solution with A = 8 have the same biases as those 
described above, suggesting that our results also do not depend 
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critically on the detailed disc wind model. The only other pa- 
rameter of significant influence on the predicted V s hif t is the 
electronic density to critical density ratio (N e /N a - = Nnx e /N a 
with 7V H , the total density) for a given line. If N e /N cr > 1 in 
the central parts of the jet, either because the accretion rate 
is increased (M acc oc Nh) or a tracer with low 7V cr (such as 
[S n] /16731) is used, the emissivity saturates and V s hif t becomes 
significantly lower than V^ f ° (though not as much as for MVC- 
dominated profiles). A ratio N e /N cr < 1 in the inner part of the 
jet, and thus a tracer with high 7V cr , is needed to retrieve the 
shape of V s hif t curves shown in Fig. [3] in the HVC-dominated 
case for high accretion rates. 

4. Comparison with DG Tau observations 

We compare here our predictions with observations of DG Tau 
carried out by Bacciotti et al. ( 2002 ) with HST/STIS. V s hif t were 
obtained from the MVC by cross-correlation techniques similar 
to the one described in Sect. 2. In Fig. |4] we plot observed 
Vshift (symbols) and predictions for both disc wind models with 
an equivalent beam size (0.1" x 50 kms -1 ) and at the same 
distances (deprojected) from the star. Predicted velocity shifts 
are plotted in several spectral lines, and are seen to differ by 
less than 10%. 

We first note that both the cold and warm solutions repro- 
duce well the trend of increased measured V s hif t with transverse 
radius in DG Tau. Based on the above study, we interpret this 
trend as due to projection effects only. Thus, measurements at 
^proj = 10, 20 AU < ^jet(Z) do not yield direct measures of V^, 
but only lower limits. Only data points at R pi0 j = 30 AU, close 
to the outer jet radius, would be expected to give true V^. 

Second, we note that V s hif t deduced from the cold disc wind 
model ^ 30^-0 km s are about by a factor 3-6 too high com- 
pared to observed velocity shifts of 6-15 kms~'. Steady, self- 
similar cold solutions with lower toroidal velocities (lower A) 
exist, but terminate too close to the star (Z/Rq < 10; Ferreira 
I1997i . This class of solutions is therefore definitely excluded 
for the MVC in the DG Tau jet. On the other hand, the warm 
disc wind solution with A = 13 fit very well observed V s hif t 
in the DG Tau jet within instrumental uncertainties, except for 
the 30 AU datapoints in region I and IV. Far from the jet axis, 
region I suffer from a poor signal-to-noise ratio, while region 
IV may be contaminated by non-axisymmetric structures in the 
external medium and/or bowshocks, similar to those at larger 
scale (Lavalley-Fouquet et al. 2000). 

Although the difference in V s hif t between R e = 1 and 3 AU 
is less than 25 %, we note that the latter R e value is essential 
for reproducing observed centroid velocities ^ 50 km s of the 
MVC in DG Tau (^ 80 km s 1 would be found for R e = 1 AU; 
see Fig-EJ?). This value is consistent with results of the diagnos- 
tic method of disc wind launching radii proposed by Anderson 
et al. ( 2003 ). However, launching radii Ro = 0.3-1 AU deduced 
with this method from the underestimated azimuthal velocities 
at R pm j = 10,20 AU, are only lower limits. Observations with 
a spatial beam < 5 AU and high-density tracers for ionized 
gas, favouring HVC-dominated line profiles, would allow to 
derive the full transverse velocity profile (oc 1 / V/?), and also 
more constraints on A, R\ and R e . Such a resolution is currently 
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Fig. 4. Velocity shifts as a function of radius at four depro- 
jected distances from the star, 17 to 91 AU, obtained by 
cross-correlation over the MVC. Symbols: DG Tau observa- 
tions (Bacciotti et al. 2002). Curves: Predictions for the cold 
(dashed; R e = 1 AU) and warm (solid; R e — 3 AU) disc wind 
solutions in [S n] ,16731 (light grey), [O i] /16300 (dark grey) 
and [Fe n] 1.644 fim (black), computed with the estimated in- 
clination angle of DG Tau jet (45°; Pyo et al. 2003 1 and with 
the beam size of HST/STIS (0. 1" x 50 km s" 1 ). 



achievable in the UV domain with HST/STIS and in the near- 
IR domain with the new generation of adaptive optics system 
operating on 8m class telescopes such as NACO on the VLT. 
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